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We consider a class of lattice topological field theories, among which are the weak-coupling limit 
of 2d Yang-Mills theory and 3d Riemannian quantum gravity, whose dynamical variables are flat 
discrete connections with compact structure group on a cell 2-complex. In these models, it is known 
that the path integral measure is ill-defined because of a phenomenon known as 'bubble divergences'. 
In this paper, we extend recent results of the authors to the cases where these divergences cannot 
be understood in terms of cellular cohomology. We introduce in its place the relevant twisted 
cohomology, and use it to compute the divergence degree of the partition function. We also relate 
its dominant part to the Reidemeister torsion of the complex, thereby generalizing previous results 
of Barrett and Naish-Guzman. The main limitation to our approach is the presence of singularities 
in the representation variety of the fundamental group of the complex; we illustrate this issue in the 
well-known case of two-dimensional manifolds. 
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I. INTRODUCTION 

One road to the quantization of a background-independent field theory such as general relativity is the spinfoam 
formalism. In this approach, which can be thought of as the covariant formulation of loop quantum gravity 0, Q , 
the Feynman path integral is realized as a sum of amplitudes associated to oriented two-dimensional cell complexes, 
aka foams. It is expected that this non-perturbative expansion could be free from the inconsistencies of perturbative 
quantum gravity (although outdated, [J H remain good reviews; see also [||). 

Inspired by the Plebanski's observation |6| that general relativity can be interpreted as a constrained BF theory with 
a simplqj-B-field, the construction of spinfoam models for four-dimensional quantum gravity has followed a remarkable 
pattern [3, One starts from a 'flat' spinfoam model analogous the one proposed by Ponzano and Regge with 
flat gauge connections with a compact structure group G as dynamical variables; Fourier analysis on the group then 
yields a formulation of the amplitudes as discrete sums over quantum numbers labelling the edges (1-cells) and faces 
(2-cells) of cell 2-complexes; to obtain an ansatz for a quantum gravitational amplitude, one eventually imposes certain 
'simplicity' constraints a la Plebanski restricting the range of these quantum numbers. 

In (Toj . however, it was pointed out that the amplitudes of the fiat spinfoam model are not well-defined, because of 
a phenomenon coined bubble divergences, and related to the pre sence of certain sets of faces in the 2-complex forming 
closed surfaces. This interpretation was then strengthened in where the bubble divergences were related to discrete 
Bianchi identities. With the understanding that these bubble divergences need a renormalization procedure, partial 
powercounting theorems were proved in (l2l [l3j. based on a combinatorial counting of these 'bubbles'. 

In our previous paper [Til ] , we have discussed in more detail the accuracy of the intuition that the divergence degree 
of a foam T can be inferred from its combinatorial, or topological structure. Defining bubbles as 2-cycles in the 
cellular cohomology of T, we have showed that there are cases where the divergence degree of T is indeed given by its 
second cellular Betti number (its 'number of bubbles'): when G is Abelian, or when T is simply connected. But, more 
importantly, we have also explained why such ideas are deceiving in more general situations. In short, the amplitudes 
can be reduced to integrals over the space of flat G-connections on T, which is also the representation variety of its 
fundamental group in G, and the structure of this space involves both the topology of T and the non- Abelian structure 
of G in a very non-trivial way. 

Here, we provide a finer description of these 'bubble divergences', using the tools of twisted cohomology. In particular, 
when the singularities of the representation variety of 7i"i(r) can be neglected (in a sense explained in sec. IIIII) . we 
find that the divergence degree of a foam is given by the value of the second twisted Betti number on generic flat 
connections. This result allows to extract a dominant part from the amplitudes, which we relate to the Reidemeister 
torsion of T. To illustrate these ideas, we use the weak-coupling limit of 2d Yang-Mills theory - a particular case of 
the flat spinfoam model, where these ideas were explored by Witten [lj|, Goldman [l6[ and many others. 

On the mathematical side, the flat spinfoam model has proved to be connected to a number of interesting problems 
in algebraic topology. In addition to the two-dimensional case, where it provides a very efficient way to compute the 
volume of the moduli space of flat connections [IH, remarkable topological invariants have been obtained in three 
and four dimensions on its basis. On 3-manifolds, the Turaev-Viro invariant can be seen as a regularization of 
the Ponzano-Regge model by a cosmological constant. In four dimensions, the Crane- Yetter invariant [181 ] provides a 
combinatorial way to compute the signature of the manifold. Note however that in these two examples, the compact 
Lie group of the original flat spinfoam model is replaced by a quantum group. Our analysis in this paper can be 
understood as a way to properly define the flat model without quantum groups, for a general cell 2-complexes (as 
opposed to the 2-skeleton of the dual cell complex to a triangulated manifold, as in the Turaev-Viro and Crane- Yetter 
models). 

Yet another motivation to study the bubble divergences of the flat spinfoam model is the following. As observed by 
Boulatov and Ooguri [20| , its amplitudes can be interpreted as Feynman amplitudes of a certain non-local field 
theory defined over a Cartesian power of the structure group, known as group field theory. The study of group field 
theory is considered a promising approach to four-dimensional quantum gravity [2ll |22| . One virtue of this approach 
is that by summing over Feynman 'graphs', which are really cell 2-complexes, the topology of spacetime is not a priori 
fixed, but could emerge dynamically. From this perspective, bubble divergences appear as a higher categorical version 



1 A 2-form B is simple if can be written B = e A e for a certain 1-form e. 
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of the usual ultraviolet divergences, and their renormalization, attempted in [12|,ll3jj, can be properly understood in 
the usual field-theoretic sense. 

The paper is organized as follows. In sec. UH we introduce in detail the flat spinfoam model, and the relevant tools 
from twisted cohomology. We use them in sec. IHII to compute the divergence degree of a foam away from singularities, 
and illustrate how the latter can be dealt with in the two-dimensional case. In sec. IIV1 we relate the dominant part 
of the amplitudes to the Reidemeister torsion. Our conclusion follows in sec. [V] 



II. DISCRETE CONNECTIONS AND TWISTED COHOMOLOGY 



In this introductory section, we define the flat spinfoam model, review the current understanding of its 'bubble 
divergences', and introduce the twisted cohomology induced on the foams by discrete flat connections. 



A. The flat spinfoam model 

Let us present in more details the generalization of the Ponzano-Regge [jjj] or Ooguri [2(| models which we refer to 
as the 'flat spinfoam model'. In both latter cases, one starts with a triangulated (spacetime) manifold (of dimension 3 
and 4, respectively), and uses the 2-skeleton of the dual cell complex to assign it a quantum amplitude. Both from the 
perspective of loop quantum gravity and of group field theory, however, it is natural to consider amplitudes defined on 
arbitrary oriented cell 2-complexes, whether dual to manifolds or not. We call foams such complexes. In this paper, 
we will only consider foams without boundary, which would arise in the Ponzano-Regge and Ooguri models in the case 
of closed manifolds. Let T be such a closed foam. We denote Ti (i — 0,1, 2) the set of its i-cells (vertices, edges and 
faces respectively), and V = E = \Ti\, F = \T2\- We also need a structure group G, which we take as a compact 
semi-simple Lie group with Lie algebra g. 

A connection on the foam T is the assignment of an element g e of the structure group G to each edge of V. These 
elements can be thought of as holonomies or 'parallel transport' operators between vertices of the foam. This is also 
the standard way to discretize a genuine connection on a cell decomposition of a manifold. The space of connections 
on r is therefore 

A = {A={ 9e ) eeTl eG E ), (1) 

The curvature of a connection A is encoded in the holonomies along the boundary of faces. It is represented by the 
map 



H :A^G* 

(2) 

' /er 2 ' 



A^(H f (A)= J] 9 l J :e 



eedf 

where [/ : e] is the incidence number of the face / on the edge e, and Hf(A) is the 'holonomy' of the connection A 
around the face /. This provides a notion of flatness on the foam: the connection is flat i£| 

H(A) = 1. (3) 

The flat spinfoam model is then defined formally as the partition function of a system of flat G-connections on T: 



Z{T,G) = f dAj[ 6{H f (A)), 
Ja /er 2 



(4) 



2 Throughout this paper, 1 denotes the unit element of the relevant group. 
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where dA — Y\ et z Tl dge is the Haar measure on A = G E , and 5(g) is the Dirac delta on G. Obviously, the support of 
this integral is the set of flat connections 

T = H-\l). (5) 

Since the curvature map H is smooth, J 7 is a smooth manifold of A whenever the unit 1 is a regular value of H. But 
this is generally not the case, as we will see. 

The spinfoam formalism consists in rewriting the partition function of the model defined in ([¥]) with integrals as 
sums over data labelling the cells of T and coming from the representation theory of G. Obviously, it is likely that the 
product of delta ditributions in (Q| is not well-defined (and to study this issue is the goal of the present article), but 
let us just ignore this difficulty for a few paragraphs. The spinfoam way to compute Z(T, G) starts with the spectral 
decomposition of the Dirac delta over G 

<%) = dim (p) x P {g), (6) 

p 

where the sum runs over (equivalence classes of) irreducible representations p of G with characters Xp- This assigns 
representations pf to every dual face, so that 

Z(r,G)= Yl J dA II dim ^/) X Pt (Hf(A)). (7) 

Then, for a fixed set of representations, one can perform the integrals over the E copies of the group, by tensoring the 
representations and using their orthogonality. 

To see how this goes, let us assume now that T is the 2-skeleton of a cell decomposition dual to a triangulated d- 
dimensional manifold M. Then, faces of T are dual to (d — 2)-simplices, while edges of T are dual to (d— l)-simplices. 
Since the latter have d (d— 2)-simplices on their boundary, each edge of T is shared by d faces. Thus, after expanding 
the characters on matrix elements of the group elements (g e ), one has to integrate over exactly d matrix elements for 
each of them. In two dimensions, the Schur orthogonality relation implies that all representations in ([7J are equal to, 
say, p. Introducing the Euler characteristic of the surface, x(M) = F — E + V, the partition function then reads 

Z(T,G) = £dim(p)*< M >, (8) 
p 

which is obviously independent of the chosen triangulation of M. This formula was shown by Witten 15] to be the 
most efficient way to compute the symplectic volume of the moduli space of flat connections for closed orientable 
surfaces (see sec. IIII Dl) . Also, for a given group, it is easy to see whether it gives a finite answer or not (in particular, 
it is always divergent for the 2-sphere and the 2-torus). 

Things go differently in higher dimensions. Denoting TA. p the carrying space of the representation p, the integral of 
d matrix elements over G is now given by 

dg D^\g) ® ... ® D^\g) = Pi ny{ u n ®-®H Pd y (9) 

G 

Here the right hand side is the orthogonal projector on the space Inv('H Pl <g>- • ■®%p ll ) of invariant tensors (intertwiners) 
in the product of the representations pj meeting at an edge e. This projector can be expanded onto a basis of 
intertwiners i e : <£)i_ 1 'Hp f <D, so that edges are finally labelled by intertwiners between the d representations meeting 
at each of them. For instance, if M is a three-dimensional manifold and G — SU(2), the irreducible representations 
on the faces are labelled by spins (j/)/er 2 ; an d all the edges are colored by the only invariant in Hj 1 ® T-Lj 2 ® ~Hj 3 
up to normalization, the Wigner 3mj-symbol. Intertwiners are contracted with each one another among the 4 edges 
meeting at each vertex of T. This gives rise to a so called 'vertex amplitude', which by duality assigns an amplitude 
to any tetrahedron. This yields the Ponzano-Regge model @ 

z PR (r,su(2))= H(-i) 2i '( 2 jf + i)I[w™Vf), (io) 

0»/er 2 / v 
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for which the vertex amplitude is the Wigner 6j-symbol, 

w™U f ) = (-i)^ ji {j % j}. (ii) 

The Ponzano-Regge model is often interpreted as a state-sum model for Riemannian 3d quantum gravity in the 
Palatini-Cartan formulation, i.e. with possibly degenerate metrics. Since the latter is nothing but topological BF field 
theory, one expects the partition function 2p R (T, SU(2)) to be independent of the triangulation, and thus to define 
a 3-manifold invariant. However, this expectation was never made precise, because the sums over spins are usually 
divergent. 

To summarize, if the spinfoam formalism gives in two dimensions a nice way to compute which is the volume 
of the space of flat connections (JHJ) , it is more difficult to obtain analogous results in highers dimensions. It turns out 
that, in these cases, the initial expression (j4|) will be more useful to study this aspect and understand the structures 
of divergences. 



B. Foams and bubbles 



Typical bubble divergences in the flat spin foam model arise when trying to prove topological invariance, i.e. 
invariance under changes of the triangulation for a given manifold. For instance, in three dimensions, there are two 
elementary moves (together with their inverses) which enable to relate triangulations of homeomorphic manifolds: the 
2-3 Pachner move, which transforms two adjacent tetrahedra into three, and the 1-4 Pachner move, which transforms 
a single tetrahedron into four. More precisely, the 1-4 move consists in adding a vertex in the bulk of a tetrahedron, 
and relate it to the four other vertices, thus creating four new edges in the triangulation. In the dual picture, this 
process gives a 3-cell dual to the added inner vertex, whose boundary is made of four triangular faces - pictorially a 
bubble. Quite clearly, these four faces are not independent, in the sense that imposing the triviality of the holonomy 
Hf around any three of them automatically enforces triviality for the remaining face. Thus, taking the formula 
literally, this gives a divergent factor 5(1). In the spinfoam formalism, this factor arises as ^2 p (dim. p) 2 . In fact, this 
phenomenon is to be expected for all inner vertices of the triangulation. Indeed, if F v is the set of dual faces of T 
'wrapping around' an inner vertex v of the triangulation, then it is easy to convince oneself to there will be an ordering 
of F v such that 

n h 7 = ( i2 ) 

feF v 

where e/ is ±1. In short, a holonomy around a face of F v can be expressed as a product of the others. This can be 
seen as a discrete version of the Bianchi identity for the curvature of a connection. This heuristic reasoning led Freidel 
and Louapre [ll| to a simple process designed at regularizing this type of divergences. The idea is simply to remove 
the group Dirac deltas along a spanning tree of the triangulation (that is a set of edges touching every vertex without 
forming a loop). 



However, an example detailed in 23[ shows that this process may not be sufficient: the Ponzano-Regge state-sum on 
Bing's house with two rooms is not absolutely convergent even though there are no inner vertices in this triangulation 
of the 3-ball0 Furthermore, this type of divergence is specific to the 3d case: its extension to higher dimensional 
manifolds is somewhat subtler. Indeed, its natural extension in dimension d would be that every (d— 3)-simplex, being 
dual to a 3-cell, would contribute a divergent factor 6(1). But this is already false in dimension 4, as one can observe 
when applying a 1-5 Pachner move on a single 4-simplex. This consists in adding a vertex to the triangulation within 
a 4-simplex, and to link it to the five other vertices of the simplex. This generates five 3-cells, while there are only 
four redundant deltas [20j . This is because the five 3-cells are not independent: they form the boundary of a 4-cell, 
dual to the vertex added to the triangulation. This phenomenon is nothing but a discrete analog of the reducibility of 
the gauge symmetry of the topological BF field theory [24J . This interesting view of divergences on manifolds will be 
studied more precisely elsewhere. 



3 Note this example falls into the case of simply connected complexes treated in [l4| . 
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Finally, we have argued that the natural setting to define the flat spin foam model uses cell 2-complexes, where the 
notion of 3-cells (and higher-dimensional cells) is obviously meaningless. Thus, we need a definition of a bubble which 
catches the intuitive picture of a 3-ball, using only two-dimensional objects. This means looking at some relations 
between faces. Recently, a first definition of a bubble has been proposed by Gurau [25] in a purely graphical way, 
which leads him to discuss the divergence degrees of group field theory 'graphs' in terms of 'bubble homology'. We 



found it more powerful in 14|, instead of that graphical method, to use a standard algebraic notion in topology: the 



cellular homology of the 2-complex, and consequently to propose the following 

Definition 1. A bubble is a 2-cycle in the cellular homology of a foam. 

With this definition, the 'number of bubbles' is the second Betti number &2(r), that is the dimension of the second 
homology group -f^fT). Describing typical situations where the divergence degree is exactly given by this number of 
bubbles was the purpose of But it is easy to see that, unfortunately, only a limited number of situations can 
be controlled with cellular homology. Let us consider again the two-dimensional case, for the torus. Choosing the 
standard cell decomposition with two edges (the non-trivial cycles) and one face, we have 

Z(S 1 xS\G)=[ dadb 5([a,b\), (13) 
Jg 2 

where the group commutator is [a,b] = aba~ 1 b~ 1 . Clearly, the Abelian case G = U(l) is special because it trivializes 
the commutators, leading to a formal divergent factor (5(1). The simplest non-Abelian case, G = SU(2), exhibits a 
different and interesting behaviour. The condition [a, 6] = 1 constrains only two parameters among the three of each 
group elements: their class angle remains free. Indeed, writing g — exp(iip g n g ■ a) 6 SU(2), where a is the vector 
made of the three Pauli matrice^ and n g £ S 2 is the direction of the rotation, 

[exp(iip a n a ■ a) , exp^b n b ■ a)] = 1 h a = ±n b . (14) 

Thus, the three-dimensional delta over SU(2) can only be used to integrate, say, the direction of fib, which leaves a 
component of the constraint trivially satisfied. Note that this explains the divergence of ((8]) in this situation. Moreover, 
this divergence is clearly not of the type J2j(2j + l) 2 , but rather J2j(^J + 1)°- 

In comparison with the discussion at the beginning of the section, this shows that one has to carefully examine 
the linear relations which may exist between the real components of the conditions Hf(A) — 1. Also, it shows that 
divergences cannot be extracted by just looking at the cell complex: they involve some non-trivial interaction between 
the foam and the structure group. This will lead us to introduce a twisted homology, already known to be relevant for 
2d Yang-Mills theory and more recently for the Ponzano-Regge model [23[ . 



C. Gauge transformations 

We saw that a discrete connection A € A can be viewed as a collection of E group elements, A = (g e ) ee n, 
representing the parallel transport operators along the edges of T. A discrete gauge transformation h is therefore a 
set of V group elements (h v ) v& r a acting at the vertices of T according to 

h - A= (h t (e)9eh^ e) ) eeri . (15) 

When r is the 2-skeleton of the dual cell complex to a triangulation of a manifold, this is indeed the effect of a gauge 
transformation on the parallel transport operators of a genuine connection in the continuum. 

We will be interested in 'factoring out' gauge transformations to identify equivalent flat connections: those lying 
along the same orbit of the G y -action. As a first step, it is convenient to reduce gauge transformations so that they 
act at one vertex only. This can be done through the standard process of contracting a maximal tree of T, that is a 
tree touching every vertex of T whitout forming loops @ 



4 The Pauli matrices are Hermitian matrices, satisfying the commutation law [(jj, tjj] = ^i^ij k &k > f° r hj = 1> 2, 3, with e 4 ^ fc the completely 
antisymmetric tensor. 

5 This will not affect the divergence degree of the foam, which turns out to be homotopy invariant, see below. 
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Definition 2. If T is a foam with V vertices, E edges and F faces, its reduction is the deformation retract of T with 
1 vertex, E — V + 1 edges and F faces. 

In the following we will only consider reduced foams T of this kind, hence with V = 1, except when explicitly stated. 
A gauge transformation is then just the conjugation of the elements g e by a single element h 6 G: 

h- A=(hg e h- 1 ) eeri . (16) 

To evaluate the effect of a small gauge transformation, we consider the differential at h = 1 of the map 

7a : G -> G £ 



/i i— > (h g e h 11 



(17) 



It is given by 



d7A|i : ^T A G E 



' ^ (i? Se »(id-Ad fle )'. 



(18) 



eeri 



where Ad stands for the adjoint representation of G on g, and for the right translation. In matrix notation: 
Ad g v = gvg~ x and R g *v — vg. 

The kernel of djA]! is the algebra of the isotropy group £(^) °f the connection A, while its image corresponds to 
the directions along which A is changed by the group action. If O a denotes the orbit through A, then we know that 

O a ~G/((A), and T a O a = imdj A{1 . (19) 

The following cases can be distinguished: 

Definition 3. A connection A is said to be reducible if it admits a non-trivial isotropy group, or irreducible if it is 
only preserved by the center ((G) of G. 

So an irreducible connection is characterized by the fact that rf7A|i is of maximal rank, rk^^^ = dimG, and the 
orbit through it is O a ~ G/((G). For G = SU(2), the reducible connections are precisely the Abelian connections, 
such that the set (g e )eeri lives in a subgroup T ~ U(l) of SU(2), or equivalently rkc^^ij = 2. In this case, the orbit is 
the homogeneous space S 2 — G/T. Finally, trivial connections, such that A g ((G) E , are called central connections. 
They are left invariant by G itself, and the rank of c?7a|i is zero. 



D. Flat connections and twisted cohomology 



Since our amplitudes are supported only on the set of flat connections F, let us now further our insight into the 
structure of this space. One useful way to describe it is in terms of the fundamental group 7Ti(r) of the 2-complex. 
Indeed, given a cell 2-complex T, one can find a presentation of its fundamental group by retracting a spanning tree, 
as described above: then, the generators a e of this group correspond to the edges e, and there is one relation per face 
which is formally exactly the flatness condition ([3]) for the generators of 7Ti(r): 

Trx(r) = ((a e ) eeri |(n4 /:el )/er 2 = 1). (20) 

e 

Notice that the relationship between presentations of groups and 2-complexes goes both ways: a finite presentation 
of a group 7r unambiguously determines a complex T. From a single vertex, draw an edge for each generator, and 
attach the faces according to the relators@ In the following we will use this process to identify a cell-complex T with 



Note that trivial relations such as aa~ 1 = 1 must not be eliminated from the presentation of n for this duality to hold. An example 
of this issue is provided by the 'dunce hat': while (a\a?a~ 1 = 1} is obviously equivalent to (a\a = 1} as a group presentation, the 
corresponding 2-complexes, the dunce hat and the disc respectively, are not. 
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one vertex and a presentation of its fundamental group 7Ti(T). The natural question is then to understand to which 
extend the amplitude Z(T, G) does depend on the chosen presentation of 7Ti(r). We will examine this question in 
section [IIIBI and show that the dominant part of Z(T,G), defined below, is invariant under certain changes on the 
presentation corresponding to 2-deformations of the foam T. Later on, in section IIV1 we will see that the dominant 
part of i?(r, G) can be furthermore expressed, under some assumptions, in terms of the Reidemeister torsion of F, 
which is known to be an invariant of simple-homotopy. 

The previous argument shows that a flat connection on T can be seen as a homomorphism from 7Ti(T) to G, 



This space is usually called the representation variety of 7Ti(r) into G (or the space of flat G-bundles over M, when 
7Ti(r) is seen as the fundamental group of a manifold M). It is independent of the chosen presentation of 7Ti(r). In 
cases of interest, G is a classical matrix group so that J 7 is a real algebraic set. 

Since the integrand of ^ is gauge invariant, it is natural to also introduce the moduli space of flat connections on 
r, i.e. .F/G, also known as the character variety of 7Ti(r), and which is also a real algebraic set. We can distinguish 
the following situations: 

• If J-/G consists in finitely many points, then it means that the tangent space to J 7 at a given point is just the 
tangent space to the orbit of this point. For instance, this happens when 7Ti(r) is a finite group. This situation 
is very similar to some cases studied in [14j . except for some non-Abelian features. However, it can be treated 
following the same lines [lij ]. simply changing all quantites coming from the cellular cohomology with their 
equivalent in the twisted cohomology we are about to describe. 

• We will be mainly interested in the case when T/G is of dimension at least one. Generically, J- has several 
irreducible components. In the following, we will concentrate our analysis on a single component, and keep in 
mind that in the end one has to sum the contributions from the different irreducible components. So we assume 
without loss of generality that T is an algebraic variety, i.e. that it is irreducible. 

To perform the integrals defining our amplitude (jj), we also need to know the local structure of the set of flat 
connections T . In particular, the following observation is key to our analysis: T = i? _1 (l) decomposes into a smooth 
submanifod J-q of the same dimension as J- (the 'generic' connections), and a set of singular connections, which is 
of smaller dimension. The set !Fq can be identified locally using the differential of the curvature map H with the 
following 

Definition 4. A flat connection <p is non-singular if dimkerdH^ = dim J 7 . Otherwise, it is singular, and we have 
dimkercLff^ > dimT^J 7 . 

The local structure of the space F was described by Goldman [l6[ , using the notion of the Zariski tangent space. 
For what concerns us, it is enough to recall that if <p is non-singular, then 



as one would expect from differential geometry. Let us emphasize however that (|22[) does not hold if <f> is singular. 
Since the starting point of our analysis is precisely this relation, the presence of singularities in T is the main limitation 
to our approach. This said, the behaviour of dH^ in the neighbourhood of a singular point is an open question for an 
arbitrary group 7Ti(r). When 7Ti(r) is the fundamental group of a closed orientable surface, which is the most studied 
situation, the variety of representations into G has singularities, but it turns out that they do not affect the partition 
function. In any similar cases, our method produces the exact divergence degree, and we plan to present in the future 
classes of manifolds for which that the singularities of the representation variety play no role. 

Let us now introduce the cohomological language which will allow us to compute the divergence degree of a foam, 
and eventually relate its dominant part to the Reidemeister torsion. First, let us switch notation and set, for any flat 
connection 




(21) 



T^Tq = kcrdff^, 



(22) 



(T, G) = drrm and <^(r, G) = dH^ 



(23) 
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for the differentials of H and 7^ at <f> and 1 € G respectively. Now, a gauge transformation by a group element h 
changes the holonomies around each face by conjugation, 

h ■ H f (A) = hH f (A)h-\ (24) 

and therefore maps a flat connection to another flat connection. Locally, this means that the directions of the orbit 
lie in the kernel of <5i, and therefore that 

Sl(T,G) o <5°(r,G) = 0, (25) 

which can be explicitly checked. In other words, the study of flat connections boils down to the cochain complex 
C;(r,G) defined by 

<— c$(r, g) < * ;(r,G) c*(r, G) < ^ (r,G) c°(r, g) <— o, (26) 

where the groups CS, Gl and Gi are respectively g, T^G E and 0^0 Thus, the tangent to the orbit 0<f, at is the 
space of coboundaries B\(T, G), the tangent to set of non-singular connections To is the space of cocyles Zi(T,G), 
and the tangent space to the moduli space of non-singular flat connections Tq/G is H^(T,G). To summarize: 

T^To = Z^(F, G), 

T O = Si(r,G), (27) 
T^To/G) =Hl(T,G). 

For notational simplicity, we will often drop hereafter the dependence on F and G in the twisted cohomology. 

Note that, because A is a Riemannian manifold, all three cochain groups are naturally equipped with inner products. 
Moreover, it is easy to check that the Euler characteristic \4> 01 the twisted cochain complex is actually independent 
of <fi, given by 

X* = bl-bl + bl=(dimG) X (T), (28) 
where x(r) is the Euler characteristic of the cellular homology of T. 

We have seen in the section HlBl that divergences are to be expected when the faces of T are not independent (in the 
sense of the cellular homology), but also more generally when the components of the curvature application H{4>) € g F 
are not independent. Locally, this happens when H is not submersive at </>, i.e. when the rank of 6^ is not (dimG) F. 
We thus introduce the following terminology, standard in differential geometry: 

Definition 5. A connection A is regular if H is submersive at A. Otherwise it is critical. 

In particular, regular flat connections on L are non-singular. In our cocomplex, the group iJ| counts the 2-cochains 
which are not in the image of 6^: — C^/B^. A flat connection is thus critical as soon as its second Betti number 
for the twisted cohomology described above, defined by 

b\ ee dim E\ = (dim G)F - rk 5\, (29) 

is non-zero. Barrett and Naish-Guzman [23| have shown that the insertion of a certain tree-like observable in the 
Ponzano-Regge partition function makes it a well defined distribution on a patch of J-"o where all flat connections 
are regular. In general, however, we expect the partition function to be controlled by flat non-singular, but critical 
connections, and the divergence degree to be controlled by according to the intuition that divergences come from 
linearly dependent components of dH^. 

Notice that, by definition, the function <f> >— > b \ is in fact constant on Tq and equal to its minimum value. We will 
use the following notation: 



More precisely, the homomorphism (f> turns the Lie algebra g into a Z[7ri(r)]-module denoted gAd^- The complex is the complex 
®7ri(r) C*(T| 2), with coefficients in gAd^,- Its cohomology is that of V with coefficients in 0Ad - 
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Definition 6. We denote 6g = min^gjr&j, the constant value of the second twisted Betti number on the space of 
non-singular flat connections !Fq. 

We will show in the next section that (unless singularities spoil the result) this number is indeed the divergence 
degree of a foam, thereby extending the result of [23j to include critical connections as well as regular ones. Examples 
where the value of bg(T) is easily computed will be detailed in the coming sections. Let us mention here the simplest 
case: when <f> is the trivial connection, one can check immediately that the twisted cohomology reduces to the standard 
cellular cohomology with coefficients in g: Cj(T, G) — C*(T, 



III. DIVERGENCE DEGREE AWAY FROM SINGULARITIES 



In this section, we compute the divergence degree of the partition function Z(T,G) away from singularities in a 
regularization of 2(T, G) using the heat kernel on G. We also study its transformation under Tietze moves, and 
illustrate our result in the well-known case where V has the topology of a closed orientable surface. 



A. Heat kernel regularization 



The heat kernel on the compact Lie group G is the fundamental solution of the heat equation 

(d T - A) K T {g) = 0, (30) 

in which A is the Casimir-Laplace operator on G. The heat kernel is a central function. Thanks to the Peter- Weyl 
theorem, it can be decomposed over the characters of G, which form an orthonormal basis of eigenfunctions of A, with 
eigenvalues C(p): 

K r (g) = Y / (dimp)e-^ Xp (g). (31) 
p 

When t goes to zero, this goes to ^2 p (dim p)Xp(ff) = 5(g) which is indeed the usual expansion of 5(g). 

For small times r, the heat kernel is localized around zero, and its behaviour is close to the Euclidean kernel. In a 
neighbourhood of the identity, 

K T (g) ~ A? imG e-^, (32) 

where \g\ is the Riemannian distance from the identity to g 7 and 

A T = (47TT)- 1 / 2 . (33) 
We thus define the regularized partition function as 

Z T (Y,G) = I dA J] K T (H f (A)), (34) 
Ja fer 2 

The integrand is obviously gauge invariant, thanks to the centrality of K T . Within this regularization, we can define 
the degree of divergence in the limit r — > 0: 

Definition 7. The divergence degree of a foam T is, when it exists, the number f2(r, G) such that the limit 

Z'(T,G) = lim A^ (RG) Z r (r,G) 

T— >0 

is finite and non- vanishing. In this case, we call this limit the dominant part of the partition function. 



It this sense, our previous results in are a particularly simple instance of those presented here. 
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B. Transformations under changes of presentation of the fundamental group 

As explained in the beginning of section Til D[ the amplitude Z(T,G) can be seen as a function of a finite group 
presentation. From this perspective, it is a natural question to ask whether different presentations of the same group 
yield the same number. It is known that finite presentations of isomorphic groups are related by a finite sequence of 
two types of elementary transformations, the Tietze moves: 

1. The Tietze moves of type 1 consist in the addition of a new generator together with a new relation expressing 
it as a word in the old generators, and the inverse operation. 

2. The Tietze moves of type 2 consist in the addition of a new relation which is implied by the old ones, and the 
inverse operation. 

Let us see how the regularized amplitudes Z T (T, G) transforms under them. 

On the foam T, the first Tietze move means that we add an edge e*, together with a group element <? e », and a face 
/* carrying the relation <? e . = w(gx, . . . , <?_e). Quite clearly, since e* only appears on the boundary of /*, one can use 
the translation invariance of the Haar measure to show that the regularized amplitude Z T (Y, G) is invariant under this 
operation: 

f E dg e *dAK T (g e ,w(g 1 ,...,g E )- 1 ) JJ K T (H f (A)) = f dg K T (g) [dA J[ K T (H f (A)), (35) 
JGxG /er 2 /er 2 

= Z T {T,G). (36) 

In the last line, we have used the fact that the integral of the heat kernel is normalized to 1 for any r. The process of 
adding to T an edge and a face as described above is called an elementary 2- expansion of T, while the inverse move 
is an elementary 2-collapse. A finite sequence of elementary 2-expansions and 2-collapses forms a 2- deformation of F, 
and the above argument shows that Z T (T, G), hence (when it exists) the dominant part of the partition function, is 
invariant under such 2-deformations, which are particular cases of simple-homotopy equivalences. 

The case of Tietze moves of type 2 is subtler. Consider the first of them, which adds a new, redundant, relation. At 
the level of the foam, this corresponds to the addition of a new face /*, such that the holonomy around /*, Hf*(A), 
is automatically trivial when the other flatness conditions hold. An obvious formal manipulation on the unregularized 
amplitude (|4]) suggests that this should change the divergence degree of the foam, but not the dominant part of Z(T, G): 

I dA 6(H f .(A)) J] KHf{A)) = 5(1) Z(T,G). (37) 
Ja fer 2 

For this reason, it has been claimed that Z(T,G) should be a function of 7Ti(r) and G only [l^. But, when applied 
to the regularized amplitude, this transformation yields 

/ dAK T (H f .(A)) J] K T (H f (A)), (38) 
Ja /er 2 

which is not related to Z T (T,G) in any simple way. The reason is that, although the set of flat connections T C A 
is unchanged, the rate at which the integrand flows away from J- is modified. Roughly speaking, the presence of an 
additional heat kernel makes the integrand more sharply peaked on T: it decreases the typical width around T of 
the connections contributing to the integral. It follows from this that the dominant part of Z(T, G) is changed by 
Tietze moves of type 2, and does depend on the actual presentation of t:i(T). Topologically, the Tietze moves of type 
2 generate homotopy equivalences of 2-complexes which are not simple. Like the Reidemeister torsion, to which it 
is related (sec. IIVI) . the dominant part Z'(T, G) may distinguish foams which have different simple-homotopy types, 
although they are homotopy equivalent. 
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C. Divergence degree on the space of non-singular fiat connections 

Our strategy to deal with the regularized integral (|34l) will be the same as the one employed by Forman [26| to 
deal with two-dimensional Yang-Mills theory: we will split the integral over A into an integral over the space of flat 
connections T and another over the normal space to T . However, like him0 we will assume that the singularities of 
T do not contribute to the integral, and thus consider its non-singular, smooth subset J-"o only. Although we do not 
understand the scope of this assumption in full generality, we will show in the section 1111 Dl how it can be checked in 
the two-dimensional case. A counter-example in three dimensions is given in the Appendix. 

On the manifold of non-singular connections J-"o, the following property holds. For each </> G J-o, the tangent space 
splits as 

T^A = T^To © N^Fq, (39) 

where T^Fq — ker<5^ and N^F^ = (ker dH^)- 1 . This property is the relevant condition to apply a generalized Laplace 
approximation to the integral in (j3"4")l . Physically, T^Tq is a local version of the space of solutions to the flatness 
equation of motion, while the normal space N^F represents the degrees of freedom of a connection which are fixed by 
the flatness condition, or equivalently which need to be fixed in order to localize the integral on flat connections. As 
far as divergences are concerned, the intuitive idea is that divergences are likely to occur as soon as the dimension of 
N^T (i.e. the number of degrees of freedom) is smaller than the number of constraints imposing flatness. 

Since the heat kernel on G decays exponentially away from the unit, we can begin by restricting the integral (134[) 
to a tubular neighborhood of JFq, equipped with the normal fibration induced by (|39[) . Assuming this neighborhood 
is small enough, we can then use the exponential mapping to pull the integral in the direction normal to J-"o to the 
normal spaces N^JFq, and consequently write H(A) = .ff (0, exp^y)), with y £ N^J-q. Denoting vo1jf the volume 
form on J-~q induced by the Riemannian metric of G E , and vo1a^.f the volume form on the normal fibers, the Haar 
measure on A can be written 

dA = vohr- voLv,^ . (40) 

From this argument, it follows that 

Z T (T,G)= f dA l[K T (H f (A)), (41) 
Ja fer 2 

/ vol,, / v„l AWl exp(- £ g^pMg Y (42) 



T-tO 



/er 2 



In order to take advantage of the Gaussian behaviour of the heat kernel K T at small times r, the Riemannian distance 
between Hf(<p, exp^(y)) and the identity can be expanded around the flat configuration (f>. This gives 

l^^exp^y))! 2 = \\dH m (y)\\l + 0(y 3 ), (43) 

It should be emphasized that the differential of Hf is only contracted with the variables y of the normal space at (p. 

We can then perform the Gaussian integral over the fibers, choosing an arbitrary basis d^ = (d 1 ^ a ) a of N^Fq, and 
write y = y a d^ a so as to integrate over the variables y a . This results in the following expression, whose terms shall 
be explained hereafter, 



/ \\dH^y)\\ 2 F \ dim{kcldH ^ voljW^) 



More exactly, since the exact value of the partition function was known independently thanks to Witten's formula Jg}, Forman knew 
that the singularities did not contribute. In our more general setting, unfortunately, no such result is available. 
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First, vol/v^o (^J>) denotes the volume spanned by the vectors of the basis in the tangent space N^F®. It comes 
from evaluating the volume form voLv^o at the saddle point y = in the basis d^. Second, as the Gaussian integral 
is performed with respect to the variables y a , the square root of the determinant of the Hessian ||dfl'^(j/)||^ in the 
basis d\ appears. In geometric terms, this is the volume spanned by the vectors (dH^(d^a)) a in the image B\. Third, 
notice that the right hand side of (j44|) is naturally independent of the choice of dj,. We give some further details on 
how to compute these quantities in practice in the next section IIII Dl 

Coming back to the full expression of Z T (T,G), the dependence on r is now extracted from the integrals, and we 
are left with 

Z r (T,G) ~ A?W 7'^> vol,.. ,45) 

Collecting the exponents of r coming both from the asymptotics of the heat kernel and from the Gaussian integral, the 
degree of divergence is G) = (dimG)F — dim(ker dH^)- 1 , with the right hand side computed on non-singular flat 
connections. To make the link with the twisted cohomology previously introduced, recall that dH^ is the coboundary 
operator 6^, hence the dimension of (ker dH^,)- 1 - is the rank of 6^. Since ((dim G)F — rk SYj = 6§, we have proved that 

Theorem 1. Whenever the ratio in (|45l) is integrable with respect to the Riemannian volume form on the set of 
non-singular fiat connections Fq, the divergence degree of a closed foam is given by 

n(r, G) = &g(r,G), 

where 6§(I\ G) is the value of the second twisted Betti number on Tq. 

This is the main result of this paper, and can be read as a precise realization of the idea that divergences may occur 
when the rank of the linearized system of flatness constraints is smaller than the number of constraints. Let us now 
illustrate this result in the two-dimensional case. 



D. SU(2) Yang-Mills theory on closed orientable surfaces 

The case when T is the cell decomposition of an orientable closed surface is well-known: it is the weak-coupling 
(or small volume) limit of two dimensional Yang-Mills theory [26|. Since the seminal works of Atiyah-Bott [27| . 
Goldman [l6[ and Witten .15] , the structure of the moduli space of flat connections has been comprehensively studied, 
and elucidated. Here, we illustrate how to use the local twisted cohomology to extract divergence degrees. For 
concreteness, we set G = SU(2) in this section. We show that there are no divergence for surfaces of genus greater 
than 2. In addition we treat the special case of the 2-torus, where non-singular flat connections are critical: we 
compute its divergence degree away from singularities, and explain why the latter do not affect the powercounting 
result. 

Theorem 2. If the closed foam T g is the standard cell decomposition of a closed orientable surface E 9 of genus g 
with 2g edges and one face (see below) , its global degree of divergence is given by 

h if .9 = 

Q(T g , SU(2)) = b&T g , SU(2)) = I 1 if g = 1 

[o if.g>2. 

The simply connected case, g = 0, is somewhat easier. We have shown indeed in [l4| that for any simply connected 
cell 2-complex T, the degree of divergence is given by f2(r, G) = (dim G) b 2 (T), where b 2 (T) is the second Betti number 
in the standard cellular cohomology of T. Moreover, it is already known for g > 2 that lim T _ ! .o Z T (T g , SU(2)) is finite, 
and is given by the integral of the combinatorial Reidemeister torsion over the moduli space Jo/SU(2) (We will 
show in the next section that this is still true for a generic foam - again, under the assumption that the singularities 
of the torsion are integrable.) 
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Let us consider the cases g > 1 in more detail. As well-known, a cell decomposition of an orientable close 2d surface 
of genus g can be reduced to a flower graph with 2g edges supporting only one face. This corresponds to the following 
presentation of the fundamental group: 

i"i(r 9 ) = (a u bx, ...,a g ,b g | [aiA] ■ • • i a g, b g] = !)• (46) 

It has 2g generators, and one relation. Here the square brackets denote the group commutator in 7Ti(T g ): [a, 6] = 
aba~ 1 b~ 1 . The space of discrete SU(2)-connections is simply obtained by mapping the generators to elements of SU(2), 
while the relator corresponds to the flatness condition: 

T g = {(a 1 ,b 1 ,...,a g ,b g ) e SU(2) 2 ^, [a u b±] ■ ■ ■ [a g , b g ] = l}. (47) 

In this two-dimensional case, it is not difficult to convince oneself that there are singular connections in T g , as the 
differential of the flatness relation is clearly not of constant rank. This means that the assumptions used in section 
IIII CI fail to be satisfied globally on F g . Also, different isotropy groups for the action of gauge transformations are 
involved, making the moduli space T g j SU(2) ill-defined as a manifold. But of course, the reasoning of the section 
IIII CI still applies on the smooth set of non-singular representations of 7ri(r s ) into SU(2). 

Actually, this two-dimensional situation admits a specific structure which simplifies the analysis: the non-singular 
flat connections are exactly the irreducible ones. We know that a G-manifold is stratified by the action of G (and 
admits an open and dense principal stratum). Here, it means that SU(2) 2s decomposes into a finite number of smooth 
manifolds, the strata, each of them consisting in the set of points whose isotropy group is conjugated to a given 
subgroup H of SU(2). Then, the simplification comes from Poincare duality which gives in 2d: 

Hl~ H l ( 48 ) 

for any flat connection <f> [l6| . This means that we have a stratification of J- g according to the different possible 
isotropy groups, coinciding with a stratification of T g according to the rank of the coboundary operator 5^ (since 
rkSh = (dimG) — bV). In particular, the irreducible flat connections, when they exist, are regular (and thus non- 
singular), whereas reducible flat connections are critical. 

To be more explicit, let us distinguish the cases g > 2 and g — 1, the latter being somehow more singular. 

• g > 2 - We can write: 

^>2=^UJ(T)UJ (G) , (49) 

where J-"o is the principal stratum: 

Jo - G T, C(4>) = C(G)} = {</» G J 7 , vkSl = dimG - 3}. (50) 

Its stabilizer is the center £(G) of G and the orbits are isomorphic to G/C(G). The space F(h) f° r H = T,G are 
smooth manifolds for which the isotropy group H is either the torus T — U(l) or G itself: 

F [H) = Q{cj>) ~ H] = {0 e J 7 , vkSl = dimG/tf }. (51) 

Now let us see how this can be used in more details. Let us first show that there are indeed non-singular 
connections where the rank of Si is 3, so that Tq is not empty. To this purpose, pick elements ai,bi G SU(2) 
trivializing the commutators [a,i, bi\ = 1 for i = 1, . . . ,g. The group elements and bi should thus have the 
same axis of rotation hi £ S 2 , but can have different axes for different i. Pulling 5^ back to the unit so that it 
acts on the Lie algebra g E = T t G E , we have: 

g 

Sl(u,v) =53(1 - Ad bi ) Ui - (1 - Ad a! ) Vi . (52) 

i=l 

for algebra elements u = (wi)j=i,..., 9 and v = (Ui)j=i,..., g - Identifying m and Vi with 3- vectors, the adjoint action 
Ad a; is just a rotation around the axis hi. Thus for each i, it is natural to decompose Ui and Vi into orthogonal 
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and parallel components to the direction hf. Ui = Ujii + Uix and similarly for Uj. Then all parallel components 
disappear, since (1 — Ad ai )vi = (1 — Ad ai )vi±. Moreover, the latter quantity still belongs to the orthogonal 
plane to Hi, where (1 — Ad) is invertible. This means that, varying Ui and Uj, each term of the sum in the above 
equation spans the orthogonal plane to n,. Finally, if at least two axes among the g directions are distinct, 
then the span of 5^(ui,Vi) is the whole algebra g. So there exist flat connections where the curvature map is 
submersive, i.e. the set of (non-singular) regular flat connections is not empty. Consequently, by the implicit 
function theorem, the latter form a smooth manifold of dimension 

dim J"o = dimker Jjr = 6g — 3, (53) 

and £>§(r g , SU(2)) = 0. Then, the reasoning of section UlI CI leads to the conclusion of a finite result away from 
singularities. We refer to the work of Sengupta [28] for details on why singularities do not contribute to the 
final result. Let us simply say that the Abelian connections correspond to taking the same axis n for all group 
elements, so that F(t) ls a manifold of dimension 2g + 2, on which rk<5^ = 2. In particular, observe that the 
relation dim T = dim ker 6^ does not apply to the singular strata. 

The most efficient way to compute the partition function in this situation is to expand the heat kernel onto 
representations and integrate the group elements using the orthogonality relation of matrix elements: 

Z r (r s ,SU(2))= f JJda.dh K r ([a 1 ,b l ]---[a g ,b g ]) = (VolG) 2g ]T +^- 2 - ( 54 ) 

It is well-defined for g > 2 when r goes to zero: 

Z T=0 (r g ,SU(2)) = (VolG) 29 ^>-( 2 ^ 2 ). (55) 

n>l 



• The torus case, g = 1 - Now we look in details at the case of the 2-torus, which is manifestly divergent if we try 
to use the above formula directly. The reason is that the isotropy group of any flat connection never reduces to 
the center of G, and is at least a U(l) torus. Through Poincare duality, this implies that the curvature map is 
not submersive on flat connections, so that non-singular flat connections are critical, i.e. H^(T g= i) ^ 0. Now, 
let us show that the the result of the section UlI CI apply nevertheless. The curvature map is here a single group 
commutator, (a, b) n- [a, b] = a6a _1 6 _1 , so that flat connections consist in rotations with the same axis h £ S 2 : 

J F i = {(«,&) € SU(2) 2 , a — exp(iip a n ■ <?), b — exp^iipt, n ■ <?), with h £ S 2 , (tpat ipb) G [0> 7r ) 2 }- (56) 

Notice the '±' in b due to the Weyl symmetry. Clearly, J 7 is a smooth manifold of dimension 4. The key point 
is that the tangent space to T is given by the kernel of (5^ whenever either a or & (say b) is not in the center of 
G. Indeed, 

Sl = d([a,b]) = (1 - Ad b ) da a' 1 - (l - Ad a ) db ir 1 . (57) 

The operators (1 — Ad a ) and (1 — Adb), seen as linear maps on q ~ K 3 , have the same one-dimensional kernel, 
the direction parallel to n: this corresponds to variations of the angles ip a , ipt, for a fixed n. Restricted to the 
orthogonal plane, (1 — Adf,) is invertible, since Ad;, is a non-trivial rotation. Thus, the equation for the kernel 
of 8\ reads 

Ua± = (l-Ad b )- 1 (l-Ad a )u b± , (58) 

and therefore fixes 2 real components of (u a ,Ub) € R 3 x R 3 (i.e. rk<5^ = 2). It expresses the condition that the 
connection remains flat under variations of the directions of a and b. We can thus evaluate the second twisted 
Betti number on the set of non-singular flat connections: since the coboundary operator 8\ is of rank 2, we have 
bl = 3 - 2 = 1. 
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But again, let us stress the existence of singularities, although T is smooth: if a, b € ((G), 5^ is the zero map. 
There it is clear that the tangent space to T is not the kernel of 5^: T^T ~ M 4 7^ ker<5^ = T$ SU(2) 2 . These 
situations correspond to an isotropy group which is G itself, and we have the stratification 

Fg=l = - F (U(1)) U -F(G) 1 (59) 



with 7( G) = C(G) 2 C ^(ua)) ■ 

Let us now apply the method described in the previous section to compute the small r behaviour of the partition 
function 

Z r (r g= i,SU(2)) = / dadb K T ([a,b}). (60) 

^SU(2) 2 

We parametrize the directions h a and n& with spherical angles (9 a , <p a ) and (#&, (^b), so that the Haar measure is 
da = sin 2 ijj a sin# a dip a d6 a d(p a , and similarly for db. The saddle points (f> G J 7 are simply given by = ±n a , and 
for now we focus on the fib = +fi a component of J 7 . To use our formula (|45[) , we need to find a basis d\ of the 
orthocomplement of ker^. First, it is convenient to change basis in the Lie algebra from the standard Cartesian 
basis ia = (ia x ,ia y ,ia z ) to the spherical basis: t„ = n a ■ id, Tg = cos9 a {costp a io x + s\mp a ia y ) — sra9 a ia z , and 
t v = — sin ifaicx + cos ip a ia y . Then, the Maurer-Cartan 1-form reads 

daa^ 1 = r n dijj a + sin-0 Q (cos^a t 8 - sin?/> a t v ) d9 a + smij) a (sin-0 a Tg + cos-0 a t v ) sin6» a d(p a , (61) 

and similarly for dbb~ l . Using this expression, we can compute the action of 5\ on tangent vectors 

y ee y°'d 0a + ^-8 Va + y e -de b + ^d Vi , (62) 
sm9 sm9 

with the vector fields evaluated on T . Here, 9 is the common value of 9 a and 9b on T . This gives 
8\(y) = 2sini/i a sini/i b |^sin(-0 a +-0 b )r e + cos(-0 a +ipb)r v ?J (y 6a - y° b ) 

+ (sin(V>„ + Vb) t v - cos(^ Q + fo) r fl ) (y^ - y^) | . (63) 

This is a rather simple operator, namely the composition of a rotation and a homothety. Its kernel is generated 
by the vectors (dg a + dg b ) and (d Va + d Vb ), in accordance with the fact that the variations of the directions fi a , 
fib must be identical in order for the connection to stay on J- . An orthonormal basis of (ker^)^ is then 

,\ _ { sin^a sin -06 / 1 1 „ \ sin^a sinV>6 ( 1 „ 1 „ 



\/sin 2 -0a + sin 2 ip b \sin tp a " sin tp b J sin^Vsin 2 ip a + sin 2 ip b \ sin V'a ° sin ip b 

(64) 

From this, we can compute the volumes entering our formula (|45l) for non-singular flat connections: 



voW^(4) = L ( 65 ) 
vol B 2 (^(4)) = 4 (sin 2 ^ a + sin 2 ^) . (66) 



As we anticipated, the determinant coming from the Gaussian integral (|66j) is singular when a and b both 
approach the center of G, where sin ip a = sin-0;, = 0. However, this singularity is cancelled by the induced 
Riemannian volume form on J- 

voljr = (sin 2 ip a + sin 2 ipb) dipadipb sin 9 d9 dip, (67) 
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Moreover, gauge transformations leave the class angles ip a ,il>b invariant, and act as rotations of the direction 
h = (8, <p). It follows that 

Z T (r ff= i,SU(2)) ~ % A T / vol S 2 / di/j a dip b = 2ttA t / #„# fe . (68) 

Note that we have multiplied the whole expression by a factor 2 to take into account the set of flat connections 
where fib — —n a instead of hb =n a , which we assumed implicitely above. 

Following Witten [l5j], this shows that the symplectic form on the moduli space of flat SU(2) connections on the 
2-torus is just dip a A dipb- Moreover, the comparison with the prediction coming from the twisted cohomology 
is successful, since we have found Sl(r g= i, SU(2)) = 1 = 6§(T s= i, SU(2)). This happens in spite of the vanishing 
of the determinant of the Gaussian form in (|66|) . which is compensated by the measure voLv .F(di) voljr. 



IV. RELATION TO REIDEMEISTER TORSION 



Our main result, (|45j) . shows a factorisation of a regulator-dependent part from an integral which is independent of 
r. Thus, our powercounting is true as soon as the integral gives a finite number. This may not be the case due to 
the singularities of J 7 , i.e. points where dimker^ > dim J". Thus, we would like to look a bit more precisely at this 
integral, and try to relate it some known quantities. We know from Witten [l5| that on a closed orientable surface, 
it is the integral of the Reidemeister torsion for the twisted cohomology. In three dimensions, Witten has shown that 
the partition function for Riemannian quantum gravity (with degenerate metrics) , of which the Ponzano-Regge model 
(IIOP is the spinfoam quantization, can be cast into the integral of the Ray-Singer torsion of the spacetime 3-manifold 
|29j . As the analytic Ray-Singer torsion is the same as the combinatorial Reidemeister torsion, we expect at the 
discrete level that the integral over fiat connections can be reduced to the integral of the Reidemeister torsion. This 
is what has been achieved by Barrett and Naish-Guzman in [23j | when the Ponzano-Regge amplitude is finite. Let us 
also mention the work of Gegenberg and Kunstatter [30j . which shows the relation between the partition function of 
the topological BF field theory in any dimension with the torsion of the spacetime manifold. Here, we reach a similar 
conclusion, but as we do not have spacetime manifolds, we find that the relevant invariant is the Reidemeister torsion 
of the cell 2-complex T, which is indeed an invariant of simple-homotopy of T. To prove this result, we will assume 
that all non-singular connections have the same isotropy type. 

As for the problem of the singularities of T and the possibility that the integral be infinite, this is related to the 
difficulty of integrating the torsion in the generic case. This problem is known in mathematics, and at the present day 



only limited results have been obtained 31] 



A. Extracting the torsion 

We first introduce a convenient and standard notation. If b and b' are two bases of a given vector space related by 
the matrix M: b' a = M a ^bp, we denote the determinant of M by [b'/b]. 

Pick up a basis c| of — g F . Then we identify the space with the orthocomplement of the image of in g F : 

Hl^^/Bl^^, (69) 

and choose a basis h 2 ^ of (B^ > ) ± . From h^ and the image of the basis d^ via 5^, we obtain another basis of written 
S^(d^) h^. The change of basis is given by 

r| = [Sl(dl)hl/4], (70) 

Thus, to compute the volume spanned by 5^, we can first compute the volume spanned by c| in g F , change basis to 
S^id^h 2 ^ and then divide by the volume of fv\: 

1 _ 1 vol (Bg)4^) 

-ol B M(dD) rj voV( c p ' ^ > 
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If <j) is regular, i.e. if — 0, then there is no hi, and vol(£2)j_(/i|) should be replaced with 1 in the above formula 
(see HI). 

Let us now describe more carefully the space of flat connections, with the motive of integrating the orbits of the 
group action. We assume that the isotropy groups of all non-singular flat connections are isomorphic, turning J-o into 
a fiber bundle. The isotropy group £(</>) of <j> is generated by keri5?. The orbit through <fi is isomorphic to G /(((/>) 
and its tangent space is the image of 5^: T^O^ = B^. We now rewrite the volume form on J-q in order to make the 
splitting between the orbit directions and their orthogonal directions in T^J-q explicit. We pick on the one hand a 
basis cj, of — T r j > G E , and on the other hand we complete the basis d 1 ^ to get a second basis on C^. This can done 
by considering a basis of (Z®) ± , and its pushforward #2(<fi) in B^. Let then h 1 ^ be a basis of the first cohomology 
space H\ = T^Fq/T^O^. By identifying with the orthocomplement of B^ in T^Fq, we can lift to a basis of 
(B^ ® Z^) 1 - . This way we get a basis, 5®(d^) d^, of corresponding to the decomposition 

T^^Bjffi^ffiZj) 1 ®^, (72) 

which can be compared to cj, through the determinant 

rl = K(4)hld),/cl]. (73) 

Furthermore, let S^(d^)* and (h^)* be the dual bases to <5°(rf°) and h^. They induce volume forms on the orbits 
and on the moduli space, denoted /\ S®(d^)* and A(^J>)* respectively. Equipped with these bases, we can rewrite the 
volume form on as 

Here, as before, the quantity vo1j^ C jb(cJ i ) refers to the volume spanned by c\ in the tangent space T^G E . This results 
in the following expression for the dominant part of the partition function: 

r T \ vo1 t , G e(cJ,) vo\( B 2 )± (hl) 

J Fa T <j> vol 5 F \ c <p) 

To simplify the following discussion, we will distinguish the reducible and irreducible cases, although the latter can 
be read as a particular case of the former. 

B. Irreducible connections 

In this case, the stabilizer is the center of the group, £(</>) ~ C(G), and via the gauge transformation map defined 
in (jTTJ) the orbits are isomorphic to G/£(G). Hence, we can pull back the integral over each orbit to an integral 
over G/((G). Since the integrand in (|75)l is gauge invariant, it is sufficient to evaluate the Jacobian of the change 
of variables at the unit in G/£(G). There, the kernel Z® of the linearized gauge transformation operator <5" is zero, 
and d°± forms a basis of the full Lie algebra G" = g. Moreover, the volume form on the orbits /\ <5?(d^)* is the push 
forward of the volume form f\(d^)* induced on g, which itself can be related to the Riemannian volume form vol B /^(G) 
according to 

A (J0\* _ vol G/C(G) _ 1 vol G/c(G) , , 

where in the last equality we introduced a basis c° of q and the corresponding determinant = [d^/c^]. Hence, 

^ (r,G) = ^ o/{(o) v„i G/c(e) j -j .g-^ A^*) ■ cm 
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At this point we can identify the Reidemeister torsion of the cochain complex G^ 

t0I > = 70^2 > ( 78 ) 

which by construction depends on the the bases of G^ and h\ of H^, but not on the bases dk. As well-known, it 
defines an invariant of simple-homotopy equivalences of the 2-complex T |32j |. 

Finally, the integral over the orbit G/((G) simply gives a factor Vol(G)/#£(G), and we obtain 

Z'(T G) V0l(G) f tor, ^'^^i) 1 ^ A (hi)* ( 79 ) 
2 (r ' G) - #C(G) i 0/G tOT ^ vol B (c«)vol gF (4) A(^) ■ (79) 

The integrand is naturally independent of any choice of basis. The last step to simplify this expression is borrowed to 
(23j . and consists in choosing the bases in a suitable way. Let I be a basis of the Lie algebra g, and choose A = I 
together with A made by one copy of I for each face. This implies vo1 8 f(c 2 ) = (vol fl (7)) F . From I, one also gets a 
basis of q e with one copy of I for each edge. Then push forward this basis isometrically to G^ = T$G E , and choose 
cj, as the result of this operation, so that voL^ge^) = (vol g (l)) E . Since the Euler characteristic of the cell-complex 
T is x(r) = F - E + 1, this gives 

2'(r,G) = (vol fl (0)" x(r) / tor voW K (^) (80) 

A particularly simple choice for I, of course, would be to pick an orthonormal basis, in which case vol g (l) = 1. Moreover, 
we are free to choose for fv^ orthonormal as well, so that vo1( B 2)_l (h 2 ^ ) = 1, and thus 

2 ' (r - G) = i?iX„ /G te *AW- <w 

This concludes our argument in the irreducible case. 



C. Reducible connections 

The reducible case proceeds along the very same lines, except for one complication: the presence of non-trivial 
cohomology classes in if 9. To obtain a basis of the Lie algebra C°, we should therefore pick a basis and via the 
identification ff9 ~ (Z9J lift it to complete the basis If c° is another basis of g, the change of basis r9 is now 
defined by 

4 = [hl4/c% (82) 

and the relation between the Riemannian volume form on G/T, where T = C(</0 is the common stabilizer of non- 
singular flat connections <j>, and the form f\(<A)* induced by cR now reads 



mY . J*ZZ_. (83 ) 

vol fl (c9) 



Trivializing locally the bundle G — > G/T, we can as before pull the integral over the orbits Q$ back to G/T: 

*M-f ([ voi e/T ) 4 AW)'- < 84 ) 

Jto/g\Jg/t J ^ vol B (c^) vo1 s f(c 2 ) '\ 
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The definition of the torsion is unchanged, tor^ = -^n, and the integral over the homogenous space G/T gives a 
trivial factor Vol(G/T) 

Z'{T,G) = Vol(G/T) (vol fl (0)" x(r) / tor vol z o(/i°) vol (B?) x(^) f\(h\)*, (85) 

where as before we used bases constructed from a fixed basis / of the Lie algebra q. Choosing moreover /i° and h 2 ^ 
orthonormal, this simplifies to 

Z'(T,G)=Vo\(G/T) f tor /\(^)*, (86) 

This last equation is the generalization of ([BTj) to the reducible case. 

To summarize the results of this section, we have showed that 

Theorem 3. If all non-singular G-connections cf> 6 on the cell 2-complex V have the same isotropy type, and if the 
Reidemeister torsion is integrable on J-" , then the dominant part of the partition function of the flat spinfoam model is 
the volume of the moduli space of non-singular flat connections Fo/G with respect to the Reidemeister torsion volume 
form, in the sense of (|8T1) and (l86l) . 



V. CONCLUSION 



With a caveat related to the singularities of the representation variety of a finitely presented goup, we have described 
the structure of the 'bubble divergences' arising in the fiat spinfoam model: their divergence degree is given by the 
second Betti number in the natural twisted cohomology. This result has allowed us to identify a dominant part in the 
amplitudes, which is be given by the volume of the character variety of the fundamental group of the foam T, with a 
volume form given by the Reidemeister torsion of V . 

One could try to use these results to define properly the Ponzano-Regge model of Riemannian gravity on 3-manifolds, 
and its higher dimensional analogues, as finite topological invariants. To this aim, and also to settle the issue of 
singularities, it will be necessary to focus attention to a specific class of foams T, such as special polyhedra. One 
could then attempt to relate the divergence degree f2(r, G), and the dominant part Z'(T, G), to the topology of the 
manifold itself. What is more, this would also permit a sharper comparison with topological quantum field theory, 
where similar divergences arise as a consequence of the continuum gauge symmetry [29j . 

From the perspective of four-dimensional quantum gravity, we hope that our results will provide a useful basis for 
the study of the divergences of 'non-flat' spinfoam models such as ,7j, 8]. Whether this hope is legitimate, the future 
will tell. 

APPENDIX: ON NON-INTEGRABLE SINGULARITIES 

Example of a non-integrable singularity 

The flat spinfoam model is defined by the choice of a finitely presented group it and a structure group G. Consider 
G = SU(2) and 

7T = (a,b,h\ [a,h] = [b,h] = 1). (87) 

It turns out that the divergence on singular connections of the Gaussian determinant in (144[) is not integrable in this 
case. In other words, (tt, SU(2)) is an example which cannot be handled by the method described in this paper. 

The set of flat connections on the corresponding foam T with three edges and two faces is determined by the relations 

J r ={(a,b,fe)eSU(2) 3 , [a,h] = [b,h] = l}. (88) 
This set has two irreducible components. 
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• If h is in the center of SU(2), i.e. h = ±1, Then, a, b can be arbitrary: 

J- irrcd = {(a,b,±l), (a,b) G SU(2) 2 } . (89) 
These are the irreducible representations of it = 7Ti(r) into SU(2). 

• If ft, is not in the center, then a, b and h have to lie in a common U(l) subgroup of SU(2). These are the Abelian 
representations. If we write a generic element g = exp(iiph ■ a), with h G S 2 the direction of the rotation and 
ip G [0, 7r) its class angle, then 

Tred = {(a, b, h) G SU(2) 3 , a = exp(±z-0 a n • a), b = exp(±iip b h • a), h = exp(iiphh ■<?)}. (90) 



Quite obviously, .Fined is of dimension 6, while J-" re d is 5-dimensional. In both cases, these are the dimensions of the 
kernel of 5^ on non-singular flat connections. It follows that 

6f rred (r,SU(2)) = 3, and 6 2 cd (I\ SU(2)) = 2. (91) 

So, we may expect the divergence to be controlled by non-singular irreducible flat connections, with Sl(T, SU(2)) = 3. 
However, this can be true only if the regular-indepedent integrals in (|45[) are finite. Here, we focus on the integral 
over the reducible representations and show that it is divergent. 

Choose as coordinates on J-" re d the variables (ip a ,ipb,iph,hh), with hh parametrized by spherical angles (9h,<Ph)- 
Geometrically, we will perform a saddle point approximation corresponding to the fact that for a fixed direction hh, 
the integrals over h a , h b are sharply peaked around in^. The variables flowing away from the space of flat connections 
are the corresponding tangent vectors: y a — y 6a dg a + y (Pa d Va and yt, = y 0b de b + y (pb d tpb . So we can write 



Z T (r,SU(2)) = / dadbdh K T ([a,h])K T ([b,h]), (92) 




x sin Va sin ip b sin ip h dip a dip b diph + Z T (T, SU(2)) irrcd , 

where a factor sin 2 9^ comes from the evaluation of (sin# a sin^f,) of the Haar measure on the {+h a = +hb = hh} 
component of and the factor 4 accounts for the other components. The action of 8\ on y a and y b is similar to that 
discussed in the 2-torus example (|63[) . The determinant of the quadratic form 1 1 5^ \ | 2 in the basis b^ = (dg a ,d Va , dg b ,d Vb ) 
can be evaluated with the same tools, yielding 

^/det(|j^|| 2 , b,/,) = 16 (sm9 h sin 2 i/j a sin 2 iph) (s'm9 h sin 2 i/j b sin 2 i/j h ). (94) 

Thus, our method produces the following tentative asymptotic equivalent for Z T (T, SU(2)) re( j: 

1 A? / j yQl \ j d^ a d^ h 
4 \J s 2 J J{o.tt] 3 sm iph 

The integral over S* 2 corresponds to the integral over the orbits of gauge transformations, which here correspond to 
rotations of the common direction hh of the three group elements, the class angles remaining fixed. In other words, the 
moduli space in this case is parametrized by the class angles ip a ,tpb and iph- However, one can see that the remaining 
integral, over the moduli space, is divergent. This means that the correct equivalent of the partition function as r goes 

to zero is not this one, and that the scaling is not given by A r ° ' (nor by any monomial in A T ) in this pathological 
case. 
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A non-integrable singularity in a Laplace integral 

Let us conclude this appendix by a very simple example showing how a non-integrable singularity can spoil the 
Laplace estimate for saddle point integrals, borrowed from 33]. Consider the numerical integal 



z T = dxdy e T . (96) 
Jr 2 

Here, the critical set is the 'cross' {x = 0} U {y — 0}, and has a singularity at (x, y) — (0, 0). A naive application of the 
Laplace approximation would give z T cx r 1 / 2 as r — > 0, since the orthogonal space to the critical set is one-dimensional. 
However, this is not the correct estimate, which turns out to be z T cx r 1 / 2 lnr. In this case, indeed, integrating along 
the 'normal fibers' yields a non-integrable singularity: formally, 




This behaviour is stricly analogous to the case where the Reidemeister torsion in not integrable in the neighborhood 
of singular connections: it signals the breakdown of the Laplace approximation, and calls for a more sophisticated 
analysis. 
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